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Two models are presented for the description of the electron screening effects 
^O that appear in laboratory nuclear reactions at astrophysical energies. The 

^vi \ two-electron screening energy of the first model agrees very well with the re- 

cent LUNA experimental result for the break-up reaction 3 He ( s He, 2p) He, 
which so far defies all available theoretical models. Moreover, multi-electron 
(^ ■ effects that enhance laboratory reactions of the CNO cycle and other ad- 

vanced nuclear burning stages, are also studied by means of the Thomas- 
Fermi model, deriving analytical formulae that establish a lower and upper 
limit for the associated screening energy. The results of the second model, 
which show a very satisfactory compatibility with the adiabatic approxima- 
tion ones, are expected to be particularly useful in future experiments for a 
more accurate determination of the CNO astrophysical factors. 
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I. INTRODUCTION 

The screening enhancement effect in laboratory nuclear reactions at astrophysical ener- 
gies has attracted a lot of attention recently, especially after the recent accomplishments of 
the LUNA collaboration at Gran Sasso • The very low energies attained for the break-up 
reaction 3 H e { 3 He, 2p) He, which is extremely important to the solar neutrino production 
|2j , revealed the real magnitude of the problem, as the screening energy obtained in that 
experiment still exceeds all available theoretical predictions. Other low energy experiments 
of the proton-proton chain 0||(past, current or planned) still need a theoretical model that 
could account for the observed enhancement. 

On the other hand the astrophysical factors for the reactions of the CNO chain have 



been obtained by performing measurements well above the Gamow peaks ( Jl0| and ref- 
erences therein) and extrapolating to lower energies without correcting for screening, thus 
committing a notable error in certain cases as it will become apparent in this work. 

Various theoretical models have been proposed so far, some of which are in conflict with 
each other (e.g. accepting J7|, or rejecting ]TTJ] the influence of the spectator nuclei) while 
others [|TJ] were applied at a time when experimental measurements were too sparse and 
inaccurate, thus their actual validity has been obscured. 



There have even been suggestions ]T3],[14| that this discrepancy between theoretical and 



experimental results is due to an overestimation of the energy losses in the experiment, 
which cause that apparent enhancement of the cross section. The most recent relevant 
experiment |IJ reports no such energy-loss deficit which means that all energy losses have 
been taken into account, and yet the observed screening energy is still higher than the 
adiabatic limit. 

Nevertheless, the prevalent belief nowadays is that a model is needed which could give 
a screening energy higher than the adiabatic limit. Although the adiabatic limit [|T5l,|T6|l is 



generally accepted, in a recent paper [T7| a simple and efficient model was proposed for the 
study of the screening effect on low-energy nuclear fusion reactions which exceeded that 
limit. In that model, the fusing atoms were considered hydrogen-like atoms whose electron 
probability density was used in Poisson's equation in order to derive the corresponding 
screened Coulomb potential energy. That way atomic excitations and deformations of the 
reaction participants could also be taken into account. The derived mean- field potentials 
were then treated semiclassically, by means of the WKB, in order to derive the screen- 
ing enhancement factor which was also shown to be compatible with the experiment. In 
that work the screened Coulomb potentials were given without details of their derivation. 
However, a detailed derivation is necessary here before two-electron configurations are stud- 
ied such as the 2 H ( 2 H, n) He reaction with a neutral projectile or the break-up reaction 
3 He ( 3 He, 2p) He. This need arises from the fact that the conventional use of a screened 
Coulomb potential is that of a Yukawa one. Nevertheless, the Yukawa one is only an ap- 
proximation (truncation) of the complete screened potential arising from the solution of 
Poisson's equation, as it will soon become apparent. Disregarding a priori higher order 
terms can possibly induce errors, especially when the experiment takes place at astrophys- 
ical energies of a few keV. 

The layout of the paper is as follows: In Sec. II the screened Coulomb potentials for 



hydrogen-like atoms are derived in a detailed fashion, which turns out to be very useful in 
Sec. Ill, where those potentials are modified in order to account for two-electron effects in 
nuclear reactions at astrophysical energies. Notably, in Sec. Ill, the screening energy for the 
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He ( 3 He } 2p) He reaction which so far remains inexplicably |18| high is reproduced to a 



very good approximation. Sec IV deals with multi-electron effects by means of the Thomas- 
Fermi model, which enables us to derive analytical formulas for the screening enhancement 
factors for reactions encountered in advanced nuclear burning stages of stellar evolution. A 
final concluding section summarizes the novelties and the results of the present work. 

II. ONE-ELECTRON SCREENING EFFECTS 

Let us consider a hydrogen-like atom with atomic number Z\. When the wave function 
of the electron is given by \l/ n ; (r, 9) then the charge density around the point-like nucleus 

is 

Pnl (r 1 6) = -e\V nl (r,6)\ 2 (1) 

Assuming spherically symmetric wavefunctions for simplicity we can solve the equation of 
Poisson for the above charge density in order to derive a screened Coulomb potential $ (r) 
around the nucleus. Note that this potential will take into account the repulsive effects of 
the point-like nucleus, by imposing the appropriate boundary conditions: 

$(oo) = , $(0) = — (2) 

r 

The second boundary condition indicates that if a positive projectile iZ^e) is in contact 
with the nucleus {Z\e) at the center of the electron cloud then there is no negative charge 
between them to reduce the Coulomb barrier. 

Let us define the screening form factor of the screened Coulomb as a function / (r) so 
that: 

$ (r) = ^/ (r) (3) 

r 

If we insert Eq.(||) and Eq. (H) into the equation of Poisson we obtain 

d 2 fni (r) 4vr 2 

-&r- = z?\*«(r)\ ( 4 ) 

which is to be solved with the boundary conditions 

f rd (oo) = finite. , f rd (0) = 1 (5) 

Equation @ with the corresponding boundary conditions constitutes a generator of 
screened Coulomb potentials which correspond to a particular excitation and ionization 
of the atomic target. According to the quantum state of the hydrogenlike atom we can use 
Eq. (f|) in order to obtain the corresponding screened Coulomb potentials. 



Actually in the above treatment there is an implicit assumption of independence between 
the nuclear and electronic degrees of freedom. This assumption can be expressed in a 
quantitative form by the formula: 

V sc (r)=V c (r) + ® e (r) (6) 

The above formula states that the screened Coulomb potential is actually the sum of the 
bare nucleus Coulomb potential V c (r) plus the electrostatic potential energy of the electron 
cloud. Therefore one can calculate $ e (r) and then add it to the bare nucleus potential 
which simplifies the calculations especially when multi-electron atoms are considered. We 
now proceed to give an alternative approach according to the above assumption, which will 
be particularly useful in our study of two-electron effects. 

Let us assume a hydrogenlike atom Z\ in its ground state, whose electron charge distri- 
bution of Eq. (JT]) can be written: 

p(r) = p (0) exp (-r/Y ) (7) 

where r = a / (2Zi). 

The charge density at the center of the electron cloud is: 



p (o) , , 

7T \CLq J 

while the electrostatic potential of the distribution is given by the solution of Poisson's 
equation: 

i_d ( 2 d$Ar)\ = _ Anp (0) exp ( _ r/ro) (9) 

r z dr \ dr J 
Upon integration we obtain : 

c 2 A ,„ N 2 ( r \ ( , , ~ r o 



Mr) = Ci + -f -47rp(0Kexp (--) ^1 + 2-^j (10) 

The electrostatic potential $ e (r) must go to zero at infinity which gives C\ = 0. At very 
large distances r >> r , due to the spherical symmetry of the distribution , any projectile 
impinging on that cloud will actually " see" a Coulomb potential of the form: 

$ e (r>>r ) = — (11) 

r 

so that C2 = —e. Inserting the values of the parameters r and p (0) into the above equation 
we obtain the formula used without details of its derivation in Ref. HIT 

e . e /_ r \ 



$ e (r) = — + - 1 + — exp (-r/r ) (12) 

r r V 2roJ 

If a positive projectile Z 2 e impinges onto the above hydrogen-like atom the total interaction 
potential energy V (r) between the two nuclei will be due to the above electrostatic potential, 
that is Z 2 e§ e , plus the repulsive potential of the nucleus Z t e : 



^Z 2 e 2 Z 2 e 2 z 2 e 2 / r \ 

^o (r) = — — + — I 1 + — I exp (-r/r ) (13) 



where 

r? = 

2(Z 



ri = ^^FT (I*) 



The reason for replacing r with Tq is that, at astrophysical energies, the electrons move 
at higher velocities than the nuclei themselves. For example in laboratory d — D reactions 
the relative nuclear velocity equals the typical electron velocity v e = ac for E = 25keV. 
Although the above assumption is particularly valid at such low energy collisions between 
hydrogen nuclei, when reactions between heavier nuclei are considered (see next section), an 
inevitable small error is involved at intermediate energies (e.g. in the vicinity of the Gamow 
peak). As it will soon become clear the WKB treatment of the penetration factor disregards 
all effects beyond the classical turning point. Therefore, inside the tunneling region, the 
wavefunction of the electron actually corresponds to a combined nuclear molecule (Z\ + Z 2 ) 
instead of the initial Z\ atom. Of course this is an approximation and, in our model, the 
intermediate stages of the wavefunction deformations are assumed to play a minor role. 
When the electron of the target atom/ion is in an excited state, we can obtain in the 



same way the corresponding potential energy used in Ref. [17 



III. TWO-ELECTRON SCREENING EFFECTS 

Let us now assume that a hydrogen-like atom Z 2 e impinges on a resting target nucleus 
Z\e. Let us further assume that the target atom has two electrons orbiting the nucleus. In 
a Hartree-Fock approximation the total potential energy of the interaction will be: 

V sc (r) = V c (r) + V n2ei (r) + V n2B2 (r) + V nie (r) + V eie2 + V eie + V e2e (15) 

that is the sum of: a) the Coulomb potential energy V c (r) between the two bare nuclei plus 
b) the interaction between the projectile (n 2 ) and the electrons (ei >2 ) of the target nucleus, 
plus c) the interaction between the target nucleus (m) and the electron of the projectile 
(e) ,plus d) the interaction between the electrons of the target V eie2 and of course , the 
interaction (eie, e 2 e) between the electron of the projectile and those of the target. In the 
above equation only the terms associated with the nuclei will be considered functions of 
the relative internuclear distance, while the electron-electron interactions will be treated 
as perturbations which will actually raise the Coulomb barrier between the two reacting 
nuclei. We consider the following channels: 

a) The nucleus-nucleus channel 

That interaction has been thoroughly studied in most text books JL0J and needs no 
further elaboration. 

b) The atom-atom channel for hydrogen like atoms. 

In most experiments, the projectile has been considered fully stripped of its electrons 
which is the case at relatively high energies. However, when the projectile is in a neutral 
state, or at least not fully ionized, its electron cloud has to be taken into account f7|. 



For such an interaction the total potential energy can be written: 

V sc (r) = V c (r) + V n2ei (r) + V nie (r) + V eie 



(16) 



When the two electron clouds interact, their mutual ground state wavefunction must be 
antisymmetric, since the electrons are identical fermions. The spatial wavefunction is neces- 
sarily symmetric therefore antisymmetry is arranged by considering an antisymmetric spin 
singlet state. Since we still work with hydrogen-like atoms, the electron-electron spatial 
wavefunction is 



*eie (j*l, r*2) 
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(17) 



where \l/ o (~^) is actually the usual ground state wavefunction of hydrogen-like atoms. 
The electrostatic potential energy of the two electrons is: 
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where "7*1,2 are the positions of the two electrons. As it will soon become apparent the 
effective interaction between the two participants of the reaction begins upon reaching the 
classical turning point. At astrophysical energies the classical turning point is hundreds of 
times smaller than the atomic radius, therefore we can safely assume that throughout the 
Coulomb barrier the two colliding nuclei constitute a combined nuclear molecule. Under 
that assumption the above integral can be calculated |Tj| so that for Z\ = Z 2 = Z we 
obtain: 



Vr, 



hZe 2 



«0 



(19) 



That positive energy will be transferred to the relative nuclear motion, increasing the height 
of the Coulomb barrier. 

The equal charge assumption combined with the fact that during tunneling the reacting 
nuclei practically coincide with respect to the electron cloud dimensions yields: V r 



ll- 2 Ci 



V nie (r) 

On the other hand each of the two electrons is actually subject to the repulsive effect 
of a screened nucleus due to the presence of the other electron. For the combined nuclear 
molecule, we have Z t = Z\ + Z 2 , while the usual variational procedure yields an effective 
atomic number for each electron Z** = Z t — 5/16. 

Therefore for the low energy reaction of two hydrogenlike atoms in their ground state, 
with equal atomic numbers Z, the interaction potential energy is : 
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where 



2Z* 



(20) 



(21) 



In the same way we can calculate the potential energy when one (or both) atoms are in an 
excited state. 

c) The nucleus-atom channel. 

For the hydrogen-like atom target an extensive study has already been given ]I7| . How- 
ever, for a two electron target atom further elaboration is needed. In that case the potential 
energy is : 

V (r) = V c (r) + V n2ei (r) + V me2 (r) + V eie2 (22) 



If we compare Eq. ( ]TB| ) with Eq.(f22|) we see that the two potentials must be approximately 
the same. Therefore the potential energy (^) can account for the nucleus-atom channel 
with a two-electron target atom, as well. 

The penetration factor P (E) multiplied by the astrophysical factor S (E) in the s-wave 
cross section formula 

(23) 

is given by the WKB method: 



P (E) = exp 



a{E) = S f P {E) 






-Edr 



(24) 



where the classical turning point is given by 

V sc (r c ) = E (25) 

At astrophysical energies the potential energy is found to be shifted by a constant screening 
energy U e which is added to the relative energy of the collision. 

For the nucleus-nucleus channel the calculation is trivial leading to P (E) = exp (— 27rn) 
where n is the Sommerfeld parameter and U e = 0. 

For the nucleus-hydrogenlike atom channel the screening energy has already been cal- 
culated |17 | 

U e = - (Z, + Z 2 ) Z 2 - (26) 

For two hydrogenlike atoms with equal charges [Ze) we have performed numerical inte- 
grations of Eq. (^) and numerical solutions of Eq. (|25|) . At astrophysical energies where 
screening becomes important the screened coulomb potential can be safely replaced by the 
quantity 

Z 2 e 2 
V sc (r) = + U e (27) 



where: 



e 2 5 e 2 

U e = -2Z (2Z - 5/16) — + - (2Z - 5/16) — (28) 

a 8 a 



which is also the screening energy for the collision of a bare nucleus (Ze) with a two-electron 
target atom. 

For the astrophysical reaction 3 H e { 3 He, 2p) He we have Z\ = Z2 = Z = 2 and the 
corresponding screening energy obtained through the above model is U e = —338 eV. Our 
result is very close to the experimental result of the LUNA collaboration U^ x = —294 ± 
47 e V. The small difference could be plausibly attributed to energy losses and experimental 
errors had it not been for another, as yet unidentified, uncertainty source for the screening 
energies and the associated astrophysical factors, which causes a small increase in the 
uncertainty for the solar neutrino fluxes. In fact the LUNA experimental results were first 
fitted using three different approaches: a) fixing the screening energy at the value 



given by the adiabatic limit, b) allowing all three parameters (S (0) , S' (0) , S" (0)J of the 
astrophysical factor and the screening energy U e to vary simultaneously, c) using higher 
energy data to fix the parameters of the astrophysical factor while varying the screening 
energy. The two model independent methods (b,c) gave considerably different screening 
energies in that work, that is U^ b ' = —323 ± 51 eV and U^ = —432 ± 29 respectively. In 
their final article |J only the (b) method was used yielding the above mentioned value of 
jjex _ 2Q4 _)_ 47 e y t^ j s obvious that if the (c) method had been used the screening energy 
would have been higher, and the adiabatic limit would have been considerably exceeded 
again. The two methods (b,c) gave respectively \TE\ the following zero energy astrophysical 



factors: S® (0) = 5.30 ± 0.08 and S^ (0) = 5.1 ± 0.1, that is a 3.9% difference. Regarding 
the solar neutrino problem that percentage admittedly leads to a negligible neutrino flux 
uncertainty for the pp and hep neutrinos (0.1%). However the uncertainty for the 8 B and 
7 Be ones can be as high as 1.5% , which should not be disregarded. Although the author 
agrees that method (b) seems more plausible than the (c) one, a better justification is needed 
for the choice of the fitting method when such low energy experiments are considered as it 
obviously constitutes a source of uncertainty. 

On the other hand for the reaction 2 H(d,p) H , where the projectile is considered 
neutral and the target is in an atomic form we have Z — 1, therefore (Eq. (|28"1)): U e = —63 
eV. If we compare that screening energy with the corresponding one for a bare projectile 
deuteron [[i~7fl(£/ e = —54) we arrive at an anticipated result. The screening energy is indeed 
higher in the neutral projectile case but not twice as much as that of the bare nucleus case. 
That is due to the repulsive effect between the two electrons which is roughly V ee = +29 eV 
and the screening of the combined nuclear molecule by the second electron which of course 
yields an effective charge which is lower than the sum of the two separate nuclei and 
consequently a sparser charge distribution. 

IV. MULTI-ELECTRON SCREENING EFFECTS 

In the framework of the Thomas-Fermi (TF) model the screened Coulomb potential for 
a neutral multi-electron atom Z% is [^(J : 



rTF i- , Z\e 



KcM = -^W (29) 



The dimensionless function <fi (r) can be obtained by the solution of the universal differential 
equation 

d 2 <P (s) = <t> 3/2 (x) 
dx 2 \/x 



where x = r/a and the screening radius is a = 0.8853^ a . 

In our approach, instead of solving numerically the differential equation, we will make 



use of an analytic approximation of the function (j> (r) given by Tietz [^TJ. Namely: 



^(ITS? (31) 

with b = 0.536. 

At first we will consider the interaction of very light nuclei with a heavy multi-electron 
atom so that we can disregard the perturbation induced by the impinging particle to the 
average electron density of the target atom, which is the sudden limit (SL) approximation. 

If we expand the TF screened potential using Tietz's approximation we obtain the 
screened potential energy of the interaction: 

rpr? Z\Zoe 2 bZ^Zoe 2 Z-iZoh 2 e ^ / 9 \ , , 

V S T / = -±-± 2 + 3^ — r + O (r 2 ) 32 

r a a z v ' 

The third term of the above potential, as well as terms 0(r 2 ) , are negligible with respect 
to the constant screening energy shift given by the second term. For example for relative 
energies of 20 keV the tunneling region for a p + 14 N reaction begins at a classical turning 
point of r c = 500 fm. At such a distance the second term is 440eV while the third is only 
7eV. Therefore, inside the barrier, where r < r c , terms proportional to the scaled distance 
r/a practically vanish. 

Inserting the above screened Coulomb potential (|32|) with Tietz's approximation into 
the WKB integral of Eq. ( p4|) and working in the same way as with the potentials of the 
previous section we obtain the screening energy for the collision of a light bare nucleus 
{Z^e) with a neutral multi-electron atom [Z\e)\ 



3 2 

U^ = -l.2\Zl'*Z 2 — (33) 

a 



TSL _ -i 01 7 4/3^ e 



Let us now suppose that the impinging nucleus (not necessarily a light one) has been 
neutralized so that it can also be considered a TF atom. We can obtain the maximum 
screening energy transferred to the relative nuclear motion by using the formula for the 
total energy of a TF atom. In fact there are three contributions to the total energy of the 
atom: The kinetic energy of the electrons, the potential energy of their interaction with the 
nucleus and the potential energy of their mutual interaction. For a neutral TF atom with 
charge Z\e the total energy is given by pOf : 



E "^~il"-^ < 34 > 



where /i = — (0) and 



Jo \dx) 



Numerically, 



rptot 
£j tf 



-20.93 Z[ ,z eV 



(35) 



(36) 



It is therefore plausible to assume that for an adiabatic limit (AL) interaction of two neutral 
TF atoms the screening energy U^p will be the difference between the total energy of the 
combined molecule and that of the two separate atoms: 



tjAL 
U TF 



-20.93 



(Zi 



\7/3 



Zi 



7/3 



77/3' 



eV 



(37) 



16 1. For 



We can compare the results of Eqs. (gj) and (|37|) with the ones obtained in Ref. 
a reaction p +% X between a bare proton p and neutral atom %X we obtain: 

Table I. The screening energies for various proton-induced astrophysical re- 
actions as obtained through the sudden (£/f L ) and the adiabatic (U AL ) limit of 



Ref. ||16|| versus the screening energies obtained through the present TF sudden 
'Utf) and adiabatic (u£k\ limit. 



(*/#) 



U e AL 

u e SL 

U TF 

a AL 



TF 



(eV) 

(eV) 
(eV) 
(eV) 



p+ 7 3 L 
186 
134 
142 

259 



p+f 
347 
281 
281 

474 



B p + 
441 
366 
359 
592 
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C p+f 
544 
462 
441 
717 



iV p + 
653 
570 
527 
847 



is 



o 



We observe that Eq. (^) practically reproduces the sudden limit £7~ L of Ref. |]T6| , while 
Eq. (^7|) gives a higher adiabatic limit than the one [Uf L ) given in Ref. |16|| . Despite the 

fact that for the reaction 3 H e { 3 He, 2p) He the electrons involved are too few to justify use 
of the above formulas it will give a sense of the validity of the present models if we apply 
our formulas on that reaction as well. In fact Eq. Q3"5p gives U^p = 166 eV while Eq. (|3"7|) 
gives U^p = 426 eV, which are admittedly reasonable bounds. 

Regarding the acceleration effects produced by the above screening energies on the non- 
resonant wings of astrophysical nuclear reactions , we can apply the usual enhancement 



factor [22 



ASL,AL) 
JTF 



(E) = exp 



a 



(SL,AL) * 



TXTi- 



TF 



E 



(38) 



where n is the Sommerfeld parameter and E is the center-of-mass energy. We obtain for 
the sudden and the adiabatic limit respectively : 
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/#!(£) ~exp 



zl /3 zU 1/2 



IE 



72 

^ 

3/2 
(fceV) 



(39) 



and 



/#(£)~exp 



Z\7i2 \/t\ + ^2 



.7/3 



^ 



7/3' 



AV2 



3£ 



3/2 
(fceV) 



(40) 



-l . 



where A = A\A 2 (A\ + A 2 )~ L is the reduced mass number. 

The above formula will be particularly useful in laboratory experiments where such nu- 
clear reactions are involved as those of the CNO cycle and the ones encountered in the final 
stages of stellar evolution (e.g. supernova nucleosynthesis). In such reactions , it is often 
impossible to measure a nonresonant cross section in the laboratory at a sufficiently low 
energy due to technical difficulties [T(J (large background of counts, beam instability etc.). 
Therefore the associated astrophysical factor cannot be accurately extrapolated to zero en- 
ergies as required for the calculation of the effective astrophysical factor S e /f that appears 
in the thermonuclear reaction rate formulas (Ref. p3| and references therein). However, 
the recent accomplishments of the LUNA collaboration with the 3 He ( 3 He, 2p) He reaction 
inspire hope that similar low-energy experiments will soon be conducted for heavier nuclei 
(multi-electron atoms) as well, in which case the above formula will help correct the low 
energy cross section measurements. 

Let us consider for example the first member of the CNO cycle, namely the radiative 

direct capture reaction 12 C (p, 7) N. The low energy cross section is dominated by an s- 

wave resonance ( ~ ) at E cm = 424 keV, while its Gamow peak at central solar conditions is 

24.5 keV. The usually employed S e ff for that reaction is the one obtained by an experiment 

24J which measured cross sections at energies as low as E cm = 138 keV disregarding all 



screening effects. According to the present models the screening enhancement of the cross 
section at such an energy would be between 2% and 3.4% which is not a negligible correction. 
It is now obvious that in any future attempt to improve the accuracy of the extrapolation 
by lowering the energy of the experiment the proposed model of this work will be very 
useful. 

The corrections are even more important for other reaction of the CNO cycle such as 
the slowest one, which controls the energy generation of that cycle, that is U N (p, 7) O. 



Various investigations |25| (and references therein) have obtained data for that reaction 
to center-of-mass energies as low as 93 keV without correcting for screening. The error 
committed at such low energies for that particular reaction can be as high as 9% (AL). 
Finally, the corrections can be dramatic if we consider very low resonances such as the 66 
keV one of the 17 (p, a) N reaction, namely of the order of 12% (AL). Such considerable 
enhancement could presumably lead even to shifts of the resonance energies themselves, 
compared to the bare nuclei measurements. 

In figures 1 and 2, there is plotted the screening enhancement factor for the most impor- 
tant astrophysical nuclear reactions of the CNO bi-cycle with respect to the center-of-mass 
energy. Close to the Gamow peaks of those reactions (Fig.l) the predicted enhancement 
of the nonresonant cross section is already 5% (at least) while at very low energies(Fig.2) 
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2 < E cm < 5keV, the nonresonant cross section of the screened reaction is expected to be 
several times larger than that of the bare nucleus one. 

Before concluding this section we should emphasize that our results agree well with the 
ones obtained by the adiabatic assumption [|Uj according to which the screening energy of 
the reaction between atoms A and B is given by the formula: 

U* E = E el (A + B)- E el (A) - E el (B) (41) 

where E e \ (A + B) is the total binding energy of the electrons in the combined atom A + B , 
and E e i (A) , E e i (B) are the total electronic binding energies of the asymptotically separated 
atomic reaction partners A, B. For example for the l2 C (p, 7) N and the 1A N (p, 7) O 
reactions respectively, the latter model gives screening energies Uf E = 444 eV and 546 eV, 
which are bracketed by the energies derived through the present models (see table) 

V. CONCLUSIONS 

Two very efficient models are presented here which reproduce the screening enhancement 
effects that appear in laboratory nuclear reactions at astrophysical energies. The first 
model, which describes two-electron effects, relies on the Hartree-Fock approximation and 
agrees very well with the recent LUNA experimental screening energy for the reaction 
3 H e { 3 He, 2p) He, which so far remains unexplained. 

The second model is based on the Thomas- Fermi theory and yields the screening energies 
for reactions encountered in advanced nuclear burning stages of stellar evolution, where 
multi-electron effects dominate. To the author's knowledge, for multi-electron laboratory 
effects, there have never been any closed formulas such as Eqs. ([39]) and fl40|), which can 
be readily used in order to correct the cross section measurements. Moreover, the latter 
model compares well with other available theories and its use is expected to be particularly 
useful in any future attempt to improve the accuracy of the CNO astrophysical factors by 
lowering the energy of the experiment. 

Finally, the fitting method used in the determination of the astrophysical factor is 
identified here as a source of uncertainty for the solar neutrino fluxes, which needs further 
elaboration and justification. 

VI. APPENDIX 

The screened potential model approach needs some elaboration regarding its actual 
effects. In Refs. |T(J and [pl| the screening energy for a collision between the atomic target 



[Zie) and the projectile (Z 2 e) was identified as 



U. - *%* (42) 



'-(I 



where the screening radius was set equal to the radius of the innermost electrons of the 
target (R a = ao/Z-i) labeling that as the "worst case". In that model, the electron cloud is 
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assumed to be unperturbed , which is the definition for the sudden limit, and there is no 
special consideration for multi-electron effects, either. 

However, in our study of hydrogen-like atomic targets |17| we have shown that the 
screening radius R a , for the same sudden limit, is actually independent of the atomic 
number Z\ and equal to the Bohr radius R a = a^., thus obtaining a smaller screening 
energy when Z\ > 1 than the one used in Refs. |TD| and |15[| . Note that if we adopt the 



"worst case" approach for the reaction 12 C (p, 7) N, then Eq. (f42|) gives the unrealistic 
screening energy U e = 1000 eV which is beyond the present TF adiabatic limit value, too. 
On the contrary, it is easy to show our approximation is valid for most practical purposes. 
The screening term U^ that we disregarded when treating Eq. ([y|) via the WKB is actually 
proportional to the scaled relative internuclear distance r/ao '■ 

2 
U^ r) = -2Z\Z 2 — - (43) 

For a typical classical turning point r c ~ 10 _2 ao we have, throughout the barrier, r < 
10 _2 ao- Hence, in the sudden limit, the ratio of the term that we considered significant 
U e (Eq. (p6|)) to the insignificant one, given by Eq. (pf3| ) , is U e /U^> > 50Z{ 1 . Obviously 
for small Z\ , which is usually the case in astrophysical reactions, our formulas becomes 
increasingly accurate. 

ACKNOWLEDGMENTS 

This work was financially supported by the Greek State Grants Foundation (IKY) under 
contract #135/2000. It was initiated at ECT* during a nuclear physics fellowship. The 
author would like to thank the director of ECT* Prof.Malfliet for his kind hospitality and 
support. 



13 



[1] 

[2] 
[3] 

[4] 
[5 

[6 
[7 



[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[ir 

[18 
[19 

[20 
[21 
[22 
[23 
[24 
[25 



REFERENCES 

R.Bonetti et al, Phys.Rev.Lett. 82 5205(1999) 

J.N.Bahcall, 1989, Neutrino Astrophysics, Cambridge University Press 

A.Krauss, H.W. Becker, H.P.Trautvetter, C.Rolfs, Nucl.Phys.A465, 150 (1987) 

A.Krauss, H.W. Becker, H.P.Trautvetter, C.Rolfs, Nucl.Phys.A467 273 (1987) 

U.Greife, F.Gorris, M.Junker, C.Rolfs, D.Zahnow, Z.Phys.A. 351, 107(1995) 

R.E.Brown, N.Jarmie, Phys.Rev.C41, 1391(1990) 

S.Engstler, A.Krauss, K.Neldner, C.Rolfs, U.Schroeder, Phys.Lett.B. 202, 179(1998) 

T.E.Liolios, Proceedings of the Bologna 2000 Conference: Structure of the Nucleus at 

the Dawn of the Century, Bologna (2000) (in press) 

C.Rolfs, Proceedings of the Bologna 2000 Conference: Structure of the Nucleus at the 

Dawn of the Century, Bologna (2000) (in press) 

C.E.Rolfs, W.S.Rodney, "Cauldrons in the Cosmos" , The University of Chicago Press, 

(1988) 

T.D.Shoppa, M.Jenk, S.E.Koonin, K.Langanke, R.Seki, Nucl. Phys.A605 387(1996) 

G.Bencze, Nucl.Phys.A492, 459 (1989) 

J.M.Bank, L.S.Ferreira, E.Maglione, J.M.Hansteen, Phys.Rev.C.53, R18(1996) 

K.Langanke, T.D.Shoppa, C.A .Barnes, C.E.Rolfs, Phys. Lett. B369 211(1996) 

H.J.Assebaum, K.Langanke, C.Rolfs, Z.Phys.A 327, 461 (1987) 

L.Bracci, G.Fiorentini, V.S.Melezhik, G.Mezzorani, P.Quarati, Nucl. Phys. A513, 

316(1990) 

T.E.Liolios, |nucl-th/0005TJTT| Eur.Phys.J.A (2000) in press 



M.Junker et al, Phys.Rev. C57 2700(1998) 

S. Gasiorowicz, "Quantum Physics", John Wiley & Sons, Inc., (1996) 

S.Flugge, "Practical Quantum Mechanics", Springer- Verlag (1998) ISBN #3540650350 

T.Tietz, J.Chem.Physics 25, 787(1956) 

T.D.Shoppa, S.E.Koonin, K.Langanke, R.Seki, Phys.Rev.C48, 837(1990) 

T.E.Liolios, Phys.Rev.C 61, 55802(2000) 

C.Rolfs, R.E.Azuma, Nucl.Phys.A227, 292(1974) 

U.Schroeder, H.W.Becker, G.Bogaert, J.Goerres, C.Rolfs, H.P.Trautvetter, 

Nucl.Phys.A476, 240(1987) 



FIGURE CAPTIONS 

Figure 1. The screening enhancement factor fpF{E) for the most important astro- 
physical nuclear reactions of the CNO bi-cycle with respect to the center-of-mass energy 
E in the region of the Gamow peaks (Egp) as calculated for central solar conditions. 
The lower (upper) solid curve represents the enhancement of the 13 C (p, 7) N reaction 
(Egp = 24.5 keV) as calculated by the above TF sudden (adiabatic) limit. Likewise, the 
dashed curves stand for the 14 A (p, 7) O reaction (Eqp = 27.2 keV), while the dotted ones 
for the 16 (p, 7) F reaction (Eqp = 29.8 keV) . Lower and upper curves always indicate 
sudden and adiabatic limits respectively. 

Figure 2. The screening enhancement factor fpp (E) for the most important astrophys- 
ical nuclear reactions of the CNO bi-cycle with respect to the center-of-mass energy in the 
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region 2 < E < 5keV. The enhancement effect is particularly accentuated at such low 
energies. The notation is the same as in Fig.l. 
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